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ABSTRACT 

A grid of stellar evolution models for Procyon A has been calculated. These models include the best physics 
availabl e to us (including th e latest opacities and equation of state) and are based on the revised astrometric 
mass of Girard et al. (1996 ). Models were calculated with helium diffusion and with the combined effects of 
helium and heavy element diffusion. Oscillation frequencies for I = 0, 1,2 and 3 p-modes and the characteristic 
period spacing for the g-modes were calculated for these models. We find that g-modes are sensitive to model 
parameters which effect the structure of the core, such as convective core overshoot, the heavy element abundance 
and the evolutionary state (main sequence or shell hydrogen burning) of Procyon A.. The /j-modes are relatively 
insensitive to the details of the physics used to model Procyon A, and only depend on the evolutionary state of 
Procyon A. Hence, observations of p-mode frequencies on Procyon A will serve as a robust test of stellar evolution 
models. 

Subject headings: binaries: visual — stars: evolution — stars: individual (Procyon, alpha Can Min) — stars: 
oscillations 



1. INTRODUCTION 

Procyon's internal structure makes it a particularly interest- 
ing target for seismological studies. Because models of stars 
at Procyon's position in the HR-diagram have convective en- 
velopes and because we know for the Sun that the solar con- 
vection zone is responsible for driving its p-mode oscillations, 
we do expect to see p-mode oscillations on Procyon. Further- 
more, because Procyon's convective envelope is very thin, any 
g-modes that exist, possibly driven by Procyon's convective 
core, may be observable at the surface, not being completely 
damped out as they are in the case of the Sun with its more 
extensive convective envelope. 

Seismology would establish the evolutionary status of Pro- 
cyon A, which could be either in the core hydrogen burning 
phase, or in the more advanced hydrogen shell burning phase, 
after the core has been exhausted and is composed purely of he- 
lium. In the latter case, the presence of an exhausted core and a 
mean molecular weight discontinuity would leave an unmistak- 
able signature on the p-mode frequency spacings due to the ex- 
istence of mode bumping where the g-mode spectrum overlaps 
(bumps into) and perturbs the p-mode spectrum. Mode bump- 
ing affects the spacing between adjacent p-mode frequencies, 
and hence, can be discerned from non-bumped p-modes. The 
g-mode spectrum would, of course, provide even more direct 
clues as to the exact nature of Procyon's core 

The Procyon binary system consists of an F5 IV-V primary 
and a white-dwarf secondary in an 40.8 year orbit. The F5 
primary (Procyon A) is a bright, nearby star with a well de- 
termined parallax and astrometric mass. As such, it presents a 
unique target for the study of non-radial stellar oscillations. Be- 
cause of the large difference in brightness between the two stars 
in Procyon, the orbital parameters have been difficult to mea- 



sure. Strand (1951) derived a mass of 1.75Mq for Procyon 
A, a mass which was incom patible with astrophysical estim ates 
from stellar interior models ( Guenther & Demarque, 1993 ) and 
the spectroscopic \ogg. Only recently, with improved orbital 
elements and parallax, and a precise angular scale derived us- 
ing the cold chr onograph (CoCo) at the NASA IR Telescope 
Facility (IRTF) (Wang et al., 1994), has it been possible to re- 
determine the mass of Procyon A ( Girard et al., 1996 ). The 
revised mass, about 1.5 Mq , is now in agreement with the as- 

trophysical mass estimate. 

unaboyer, Demarque & Uue ntner (1VV8), stimulated by the 
SONG ground-based campaign ( Pilachowski et al., 1998 ) to de- 
tect oscillations on Procyon, made a study of the evolutionary 
status of Procyon A, using up-to-date physical input and the re- 
vised Procyon A mass. Because the SONG data analysis is not 
yet complete, only a brief summary of these Procyon models 
was published 

However, the situation has changed dramatically in the last 
few weeks, with the successful identification of individual 
mode frequencies in the giant star a UMa A by B uzasi et al. 
(1999), using the small telescope on the WIRE space mission. 
Because Procyon is on the WIRE target list for observations 
in the near future, theoretical models and their predicted fre- 
quencies have become particularly relevant. We present in this 
paper a more detailed description of our interior models, and a 
list of individual p-mode frequencies for our reference model. 
We also include a table of averaged large and small p-mode 
spacings and the characteristic g-mode period spacings for our 
models. 
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Table 1 
Model Characteristics 





Mass 


AY 




More 


M s « 


log 


<Au > 


<A^i> 


<A^2> 


<A^ 3 > 


<Su > 


<8u,> 


n 


Model 


(M ) 


AZ 


z/x em 


(M Q ) 


(M ) 


(L/L ) 


(pHz) 


(pHz) 


GuHz) 


(AtHz) 


(pHz) 


OiHz) 


0*Hz) 


Standard 


1.50 


3.03 


0.0253 


0.118 


1.17E-4 


0.85868 


53.95 


54.15 


51.38 


54.40 


4.75 


7.01 


60.58 


subgiant 


1.40 


2.94 


0.0239 


0.000 


5.62E-5 


0.85838 


51.89 


48.71 


46.23 


45.64 


6.03 


7.35 


16.87 


low L 


1 .50 


2.47 


0.0243 


0.113 


2.31E-4 


0.83747 


54.19 


54.34 


51.62 


51.80 


4.69 


7.06 


59.35 


no wind 


1.50 


3.12 


0.0196 


0.122 


1.36E-5 


0.85852 


52.99 


53.22 


50.50 


50.69 


5.09 


7.26 


62.02 


overshoot 


1.50 


2.53 


0.0248 


0.133 


1.04E-4 


0.85848 


53.80 


54.06 


51.28 


54.25 


4.64 


6.92 


64.77 


low Z 


1.50 


2.02 


0.0206 


0.107 


1.15E-4 


0.85860 


53.87 


54.05 


51.41 


54.52 


5.23 


7.10 


56.35 


high M 


1.54 


1.90 


0.0247 


0.121 


1.07E-4 


0.85854 


54.75 


54.95 


54.97 


52.32 


4.58 


7.17 


61.65 


Zdiff 


1.50 


3.01 


0.0253 


0.118 


1.19E-4 


0.85847 


53.94 


54.15 


51.38 


54.40 


4.73 


7.01 


60.70 



2. MODELING PROCYON A 

Models for Procyon A were calculated using the Yale stel- 
lar evolutio n code, in its non-rotating configuration ( Guenther 
et al ., 1992). These models included the latest OPAL opaci- 
t ies (Iglesias & Rogers, 1996), and the OPAL equation of state 
( Rogers, Swenson & Iglesias, 1996). The diffusion coefficients 
are from [Thoul, Bahcall & Loeb (19940. Most of the mod- 



els only included helium diffusion. A single model which in- 
cluded both helium and heavy element diffusion (treated as 
a single mean heavy element, Z) was evolved to study the 
effects of heavy element diffusion on Procyon A. The mod- 
els included the effects of a wind mass loss in the diffusion 
equations. A wind velocity of v w = — M-/{Airpr 2 ) was in- 
cluded in the diffusion equations, and a solar mass loss rate 
M = 2x 10 _14 Mq /yr was assumed. This wind velocity was 
large enough to effectively suppress the diffusion in the outer 
layers of the model. One model was calculated without wind 
loss. The models employed a solar calibrated mixing length, 
though the main sequence models of Procyon A have a very thin 
surface convection zone and so the structure of these models are 
largely independent of the mixing length. None of our models 
includ the effects of rotation. As Procyon A is not a rapidly ro- 
tating star, the structure effects of rotation will be much smaller 
than the other factors considered in this paper. Rotation likely 
leads to mixing in stars, an d this mixing can supre ss diffusion in 
the outer layers of the star ( Chaboyer et al., 1995 ). This effect is 
similar to the inclusion of a stellar wind (which also suppresses 
diffusion in the outer layers), and so a rough estimate of the 
importance of rotational mixing can be made by comparing the 
standard case to the no wind case. 

The fun damental properties o f Procyon A are: M = 1 .50 ± 
0.05 M Q (|Girard et al., 1996|), tt _ 0.2832" ± 0.0015 (Gi- 



rard et al., 1996), f ffi = (18.64 ±0.87) x lO^ergscnT^-1 



( Smalley & Dworetsky, 1995|) and an angula r diameter of 
(p = 5.51 ±0.05 mas ( Mozurkewich et al., 1991 ) which imply 
L = (7.22 ±O.35)L and R = (2.09 ± 0.02) Rq. The central 
values for the mass, radius and luminosity were adopted for the 
standard model of Procyon A (see below). The models for Pro- 
cyon A were evolved from the ZAMS until they reached the 
observed radius. In an iterative procedure, the helium abun- 
dance was adjusted in the ZAMS model and the model was 
re-evolved, until the model matched the observed radius and 
luminosity of Procyon A. The starting helium abundance was 
constrained such that 1.5 < AY /AZ < 4. In addition, the final 
surface abundance of Z/X was co nstrained to be within 0. 1 dex 
of the solar value Z/X = 0.0245 ( |Grevesse & Noels, 1993] ), to 
satisfy the observed constr aint that the surface abundances of 
Procyon A are near solar (Takeda, 1994; Takeda et al., 1996 



Kato et al., 1996| ; |Yushchenko & Gopka, 1996| ). Models which 
met all of the above criteria had their pulsation frequencies cal- 
culated using Guenther' s non-r adial non-adiabatic stellar pulsa- 
tion program ( Guenther, 1994| ). 

In total, eight different calibrated models of Procyon A were 
calculated and pulsed (Table [j]). The standard model (line 1 in 
Table [[]) uses the best estimate for the mass, radius and lumi- 
nosity, does not include heavy element diffusion or overshoot at 
the convective boundaries, and uses our calibrated solar value 
for the mass fraction of the heavy elements (Z = 0.018). The 
other seven models involved changing a single parameter from 
the standard model: (a) lower mass (subgiant); (b) evolution to 
a somewhat lower luminosity (low L); (c) no wind loss; (d) con- 
vective core overshoot of 0. 1 pressure scale heights (overshoot); 
(e) low heavy element abundance (low Z, with Z = 0.015); (f) 
higher mass (high M); and (g) heavy element diffusion (Z diff). 
All but the subgiant model are in the main sequence phase of 
evolution. The lower mass model is in the hydrogen shell burn- 
ing phase of evolution. 

3. PULSATION PROPERTIES 

The average first-order frequency spacing < Av > of the p- 
modes, with Av(n,£) = v(n,£) — v(n — 1,1), is likely to be the 
first quantity determined by stellar seismology. Table [j] includes 
the frequency spacing averaged over n = 5 to 25 for £ = 0, 1,2 
and £ = 3. The first-order frequency spacing is sensitive to 
the structure of the outermost layers of the star. Table [j] also 
includes the average second-order frequency spacing (< 8v >, 
averaged over n — 5 to 25) for I = and 1. The second order 
frequency spacings, 8v(n,£) = v(n,£) — v(n — l,£ + 2), are sen- 
sitive to the structure of the innermost layers of the star. Table 
[ljincludes the characteristic period spacing II of the g-modes. 

A list of modes from £ = to 3 and n = to 24 for the stan- 
dard model is given in Table [| In this table, n p is the number 
of nodes in the eigenfunction that are p-mode like (with respect 
to their phase characteristics), and n g is the number of nodes in 
the eigenfunction that are g-mode like. The gaps in the spacing s 
are due to the n-value naming convention (Unno et al., 1989). 
In regions of variable mean molecular weight, the frequencies 
of the g-modes trapped in the interior may overlap with the fre- 
quencies of the /?-modes trapped in the envelope, "bumping" 
the /?-modes, and resulting in the mixing of p-mode like modes 
and g-mode like modes in the p-mode eigenfunction. As a con- 
sequence of mode bumping the large spacing varies irregularly 
in the frequency range where mixed modes (i.e. mode bump- 
ing) occurs. In the standard model only the lowest frequency 
p-modes (lowest n) are mixed but for the subgiant model, all of 
the nonradial p-modes we calculated, up to n = 25, are mixed. 
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Fig. l.The large frequency spacing Av is plotted as a function 
of frequency v for the standard and subgiant models. Note that 
irregularities in the large frequency spacings extend to higher 
frequencies for the subgiant model as compared to the standard 
model. 

In Figure 1 we show the large spacing for the standard model 
and for the subgiant model. As just described, we see that the 
irregularities in the large spacings extend to higher frequencies 
for the subgiant model than the standard model. 

The frequencies and characteristic frequency spacings de- 
pend on the radius of the star. To first order, (5R»/R* ~ 8v/v. 
The radius of Procyon A is known to within 1 % . Thus, the 
calculated pulsation frequencies are uncertain at the 1% level 
due to the error in the radius. In addition, uncertainties in 
the modeling of the superadiabatic layer (in both the evolution 
and pulsation calculations) leads to an estimated uncertainty of 
~ ±0.5% in the p- mode frequency calculations (see Guenther 
& Demarque (1996)). The total error associated with the calcu- 
lated pulsation frequencies shown in Table [j] is approximately 
~ 1.5%. Thus, differences in Av greater than ~ 1/iHz are sig- 
nificant. From Table [l] we see that only the subgiant model has 
Av which are significantly different from the others. When the 



large spacing is observed on Procyon, it will immediately tell 
us whether or not Procyon has reached the subgiant phase of 
evolution. 

The characteristic period spacing (II) of the g-modes shows 
a much larger variation than Av (Table [j]). Once again, due to 
errors in the models and radii of Procyon A, only differences 
in II greater than ~ 1 /j,Hz are significant. In this context, it is 
clear that the detection of g-modes in Procyon A would allow 
one to (a) determine its evolutionary status; and (b) determine 
if appreciable overshoot (> 0.05 H p ) is occurring at the edge of 
the convective core in Procyon A, and/or provide an estimate of 
the interior metallicity of Procyon A. 

4. SUMMARY 

A grid of stellar evolution models for Procyon A has been 
calculated and their pulsation properties investigated. The grid 
of models was chosen to examine the effect of varying the 
physics used to construct the models (such as the inclusion of 
overshoot, heavy element diffusion or a stellar wind) as well as 
uncertainties in the observed parameters of Procyon A (mass, 
luminosity, radius and heavy element composition). The prop- 
erties of these models (including the average frequency spacing 
for p and g-modes) are listed in Table |lj while Table g presents 
a list of individual frequencies for our standard model of Pro- 
cyon A. The predicted p-mode frequencies are fairly similar for 
all models, except for the model which uses a lower (1.4Mq ) 
mass for Procyon A. This model is in the subgiant (shell hy- 
drogen burning) phase of evolution, while all other models are 
on the main sequence. If Procyon A is a main sequence star, 
the predicted frequency spacings are fairly regular, while mode 
bumping leads to irregular frequency spacings in the subgiant 
star. As the predicted frequencies are relatively insensitive to 
changes in the input physics, or the assumed mass, radius, lu- 
minosity and heavy element abundance of Procyon, the detec- 
tion of /?-modes on Procyon A will serve as a robust test of 
stellar evolution theory. Procyon A is on the WIRE target list 
for observations in the near future. As WIRE has been success- 
ful in det ecting individual mo de frequencies in the giant star a 
UMa A (Buzasi et al., 1999), it is hoped that the predictions 



presented here for the pulsation properties of Procyon A will 
soon be tested. 

This work was supported in part by an Canadian NSERC 
grant to D.B.G. 
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Table 2 

Standard Model Frequencies (in /xHz) 



£ 


n 


p 


n s 






5v 


I 






n g 




Ai> 














129.02 






2 





2 


2 


208.98 







2 


2 





227.83 






2 


1 


2 


1 


236.75 


27.77 





3 


3 





280.35 


52.52 


7.84 


2 


2 


3 


1 


272.51 


35.76 





4 


4 





330.93 


50.59 


8.96 


2 


3 


4 


1 


321.97 


49.47 





5 


5 





382.26 


51.32 


7.71 


2 


4 


5 


1 


374.55 


52.57 





6 


6 





434.86 


52.61 


6.47 


2 


5 


6 


1 


428.39 


53.84 





7 


7 





487.66 


52.80 


6.26 


2 


6 


7 


1 


481.40 


53.01 





8 


8 





541.74 


54.08 




2 


7 


7 





494.03 


12.63 





9 


9 





596.51 


54.76 




2 


8 


8 





536.95 


42.92 





10 


10 





650.31 


53.80 




2 


9 


9 





591.97 


55.01 





11 


11 





702.97 


52.66 




2 


10 


10 





646.02 


54.06 





12 


12 





755.02 


52.05 




2 


11 


11 





698.82 


52.80 





14 


13 





807.81 




4.21 


2 


12 


12 





750.94 


52.11 





15 


14 





862.40 


54.60 


4.33 


2 


13 


13 





803.60 


52.66 





16 


15 





917.64 


55.23 


4.35 


2 


14 


14 





858.07 


54.47 





17 


16 





972.79 


55.16 


4.35 


2 


15 


15 





913.28 


55.21 





18 


17 





1027.70 


54.91 


4.28 


2 


16 


16 





968.45 


55.16 





19 


18 





1081.82 


54.12 


4.21 


2 


17 


17 





1023.42 


54.97 





20 


19 





1135.90 


54.08 


4.21 


2 


18 


18 





1077.61 


54.19 





21 


20 





1 1 OA Of) 


54.39 


4. 16 


2 


19 


1 9 





1 1 3 1.70 


54.08 





22 


21 





1244.81 


54.51 


4.10 


2 


20 


20 





1186.13 


54.44 





23 


22 





1299.77 


54.97 


4.02 


2 


21 


21 





1240.71 


54.58 





24 


23 





1354.69 


54.91 


3.87 


2 


22 


22 





1295.76 


55.05 
















2 


23 


23 





1350.82 


55.06 




-2 





2 


129.55 






3 


2 


3 


1 


286.72 









1 




181.62 






3 


3 


4 


1 


339.15 


52.43 




1 


2 




236.71 


55.08 




3 


4 


5 


1 


393.34 


54.19 




2 


3 




281.78 


45.07 




3 


5 


6 


1 


447.87 


54.53 




3 


4 




302.21 


20.43 


15.49 


3 


6 


7 


1 


502.07 


54.20 




4 


5 




350.49 


48.28 


11.34 


3 


7 


8 


1 


557.30 


55.23 




5 


6 




403.82 


53.32 


10.47 


3 


9 


9 





612.62 






6 


7 




457.48 


53.66 


9.61 


3 


10 


10 





666.24 


53.62 




8 


9 




511.24 






3 


11 


11 





718.88 


52.65 




9 


10 




566.13 


54.89 




3 


12 


12 





771.10 


52.22 




10 


10 





620.65 


54.52 


8.03 


3 


13 


13 





824.58 


53.47 




11 


11 





673.85 


53.20 


7.62 


3 


14 


14 





879.65 


55.08 




12 


12 





726.14 


52.29 


7.26 


3 


15 


15 





935.03 


55.37 




13 


13 





778.19 


52.04 


7.09 


3 


16 


16 





990.31 


55.29 




14 


14 





831.69 


53.51 


7.12 


3 


17 


17 





1045.11 


54.80 




15 


15 





886.65 


54.96 


7.00 


3 


18 


18 





1099.24 


54.13 




16 


16 





941.82 


55.16 


6.79 


3 


19 


19 





1153.64 


54.40 




17 


17 





996.92 


55.10 


6.61 


3 


20 


20 





1208.28 


54.64 




18 


18 





1051.45 


54.53 


6.34 


3 


21 


21 





1263.15 


54.87 




19 


19 





1105.40 


53.95 


6.17 


3 


22 


22 





1318.42 


55.27 




20 


20 





1159.70 


54.30 


6.06 


3 


23 


23 





1373.45 


55.03 




21 


21 





1214.21 


54.51 


5.93 
















22 


22 





1269.01 


54.80 


5.86 
















23 


23 





1324.15 


55.14 


5.72 
















24 


24 





1378.99 


54.85 


5.54 















